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In lattice QCD, a confining potential for a static quark-antiquark pair can be computed with the 
Wilson loop technique. This potential, dominated by a linear potential at moderate distances, 
i—i is consistent with the confinement with a flux tube, an extended and scalar system also directly 

! observable in lattice QCD. Quantized flux tubes have also been observed in another class of con- 

^ finement, the magnetic confinement in type II superconductors. On the other hand the solution 

of Schwinger Dyson Equations, say with the Landau gauge fixing and the truncation of the series 
C\( of Feynman diagrams, already at the rainbow level for the self energy and at the ladder level for 

the Bethe Salpeter equation, provides a signal of a possible inverse quartic potential in momen- 
tum space derived from one gluon and one ghost exchange, consistent with confinement. Here 
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Q\ we address the successes, difficulties and open problems of the matching of these two different 

perspectives of confinement, the Schwinger-Dyson perspective versus the flux tube perspective. 
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1. Introduction 

Bjorken [1] asked, "how are the many disparate methods of describing hadrons which are now 
in use related to each other and to the first principles of QCD?". This is an important question, 
since the non-pertubative nature of QCD leads to different truncation schemes or models, which 
turn out to be difficult to match. Apparently chiral symmetry breaking and scalar confinement are 
conflicting, because chiral symmetry breaking requires a chiral invariant coupling to the quarks, say 
a vector coupling like in QCD. Here we show a possible solution to this old conflict of hadronic 
physics, which remained open for many years. 

This talk has two goals. Our first aim is to compute in Lattice QCD a confining flux tube, to 
illustrate scalar confinement. Because the meeting QCD-TNT is partly focused on gluons, and be- 
cause the flux tubes for mesons and for baryons have been already studied in the literature, here we 
study the novel case of the hybrid gluon-quark-antiquark system. In Bicudo et al. [2] and Cardoso 
et al. [3] it was shown that for the static gluon-quark-antiquark potential V when the segments 
between gluon-quark and gluon-antiquark are perpendicular, V is compatible with confinement re- 
alized with a pair of fundamental strings, one linking the gluon to the quark and the other linking 
the same gluon to the antiquark. However, if the segments are parallel and superposed, the total 
string tension becomes larger and agrees with Casimir Scaling measured by Bali [4]. This can be 
interpreted with a type-II superconductor analogy for the confinement in QCD, see figure 14, with 
repulsion of the fundamental strings and with the string tension of the first topological excitation 
of the string (the adjoint string) larger than the double of the fundamental string tension. Here we 
go beyond this study of the static potentials, we compute the chromoelectric and chromomagnetic 
fields and study the shape of the flux-tubes . We study two different geometries for the hybrid 
system, one with a U shape and another with L shape, figure. 1. 
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(a) U shape geometry. (b) L shape geometry. 
Figure 1: gluon-quark-antiquark geometries, U and L shapes. 

Our second aim is to show an example of a quark model for a possible coupling of a scalar 
string to quarks, that matches the apparently conflicting vector coupling of QCD with a scalar con- 
finement. The QCD Lagrangian is chiral invariant in the limit of vanishing quark masses, because 
the coupling of a gluon to a quark is vector like, and it is this vector coupling that apparently con- 
flicts with scalar confinement. We also solve the mass gap equation for the dynamical generation 
of the quark mass with the spontaneous breaking of chiral symmetry, and we indeed generate both 
the constituent quark mass and the scalar confinement. 
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2. The Wilson Loops and Chromo Fields for exotic systems 

In principle, any Wilson loop with a geometry similar to that represented in figure 2a, describ- 
ing correctly the quantum numbers of the hybrid, is appropriate, although the signal to noise ratio 
may depend on the choice of the Wilson loop. A correct Wilson loop must include an SU(3) octet 
(the gluon), an SU(3) triplet (the quark) and an SU(3) antitriplet (the antiquark), as well as the 
connection between the three links of the gluon, the quark and the antiquark. 




(a) (b) 

Figure 2: (a) Wilson loop for the gqq and equivalent position of the static antiquark, gluon, and 
quark, (b) Simple Wilson loops that make the gqq Wilson loop. 



The Wilson loop for this system is given by: 

W gm = ^Tr[X h T + X a T}Tr[X h SiX a S 2 ] (2.1) 
= W x W 2 - l -Wi (2.2) 

where W\, W 2 and W3 are the simple Wilson loops shown in figure 2b. When r\ = 0, W\ = 3 
and W 2 = W3, the operator reduces to the mesonic Wilson loop and when ji = v and r\ = r 2 = r, 
W 2 = Wl and W 3 = 3, W gm reduces to W gm {r,r,t) = \W(r,t)\ 2 - 1, that is the Wilson loop in the 
adjoint representation used to compute the potential between two static gluons. 

In order to improve the signal to noise ratio of the Wilson loop, the links are replaced by "fat 
links", 



E//I is) - P su{3) ( Up (s) + w £ U y (s) ^ (s + V) Ul (s + n) j , (2.3) 



using w = 0.2 and iterating this procedure 25 times in the spatial direction. 
We ontain the chromoelectric and chromomagnetic fields with, 



(* 2 > = ^-^> 0.5) 



and the plaquette is given by 



Ppv(s) = 1 - ^ReTr [U^s)U v (s + 11)^(1 + v)U}(s)] (2.6) 
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We only apply the smearing technique to the Wilson loop. The relationship between various pla- 
quettes and the field components is: 



component plaquette 

E x Pxt 

Ey Pyt 

E z P zt 

Bx Pyz 
By Pzx 
B z Pxy 

The energy (e) and action (y) density is given by 

e= l -((E 2 ) + {B 2 )) (2.7) 



r=\((E 2 )-(B 2 )) (2.8) 



3. Results for the fluxt tube created by a gluon, a quark and an antiquark 

Here we present the results of our simulations with 141 24 3 x 48, j8 = 6.2 configurations 
generated with the version 6 of the MILC code [5], via a combination of Cabbibo-Mariani and 
overrelaxed updates. The results are presented in lattice spacing units a = 0.072 fm. 

In the following pictures we present the results for the chromolectric and chromomagnetic 
fields and for the energy and action density. The figures 3 to 6 show the results for the U shape 
geometry and figures 7 to 10 the results for the L shape geometry. 

In figure 1 1 we present the results for the U geometry along x = 0, this corresponds to the 
axis along gluon and the middle distance between quark-antiquark, and y = 4, this corresponds 
to the axis along quark-antiquark. This results show that when the quark and the antiquark are 
superposed the results are consistent with the degenerate case of the two gluon glueball. In figure 
12 we present the results for the L geometry along the segment gluon-antiquark and along segment 
gluon-quark. And finally, in figure 13 we show a comparison between the results for (r\ = 0, r2 = 8) 
and (n = 8,r2 = 8) along segment gluon-antiquark in the L geometry, when the gluon and the 
antiquark are superposed this results are consistent with the static quark-antiquark case. 
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(a) Chromoelectric Field (b) Chromomagnetic Field 





(c) Energy Density 



(d) Action Density 



Figure 3: Results for the U geometry with d = and I = 8. In this case, the quark and the antiquark 
are superposed and the results are consistent with the degenerate case of the two gluon glueball. 
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(a) Chromoelectric Field (b) Chromomagnetic Field 





(c) Energy Density 



(d) Action Density 
Figure 5: Results for the U geometry with d = 4 and I 
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(a) Chromoelectric Field (b) Chromomagnetic Field 




(d) Action Density 




(c) Energy Density 



Figure 6: Results for the U geometry with d = 6 and I 
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(a) Chromoelectric Field (b) Chromomagnetic Field 
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(c) Energy Density 



(d) Action Density 



Figure 7: Results for the L geometry with r\ = and = 8. The gluon and the antiquark are 
superposed and the results are consistent with the static quark-antiquark case. 





(a) Chromoelectric Field (b) Chromomagnetic Field 



(c) Energy Density 




(d) Action Density 



Figure 8: Results for the L geometry with r\ — 2 and r>i — 8 
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(a) d = and / = 8 for x = (b) d = and / = 8 for y = 4 (c) d = 2 and / = 8 for x = 




-6-4-2 2 4 



(d)d = 2and/ = 8for;y = 4 




-4-2 2 4 





-6 -4 -2 



-4-2 2 4 



(e) J = 4 and / = 8 for x = (f) J = 4 and / = 8 for y = 4 (g) d = 6 and / = 8 for x = 

0.020 f 




-6-4-2 2 4 



(h)<i = 6and/ = 8for;y = 4 
Figure 1 1 : Results for the U geometry. 
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4-2 2 4 



(a) r\ = and r2 = 8, along (b) n = and r2 = 8, along (c) ri = 2 and = 8, along 
segment gluon- antiquark segment gluon-quark segment gluon- antiquark 




-4-2 2 4 



(d) ri = 2 and r2 = 8, along segment 
gluon-quark 




-4 -2 2 4 6 



-6 -4 -2 2 4 



-4 -2 2 4 6 



(e) r\ = 8 and r2 = 8, along (f) ri = 8 and r2 = 8, along (g) r\=6 and = 6, along 
segment gluon- antiquark segment gluon-quark segment gluon- antiquark 




-6 -4 -2 2 4 6 



(h) r\ = 6 and r2 = 6, along segment 
gluon-quark 



Figure 12: Results for the L geometry. 




• Action Density -> rj=8 and r2=8 
a Eneigy Density -> =8 and r 2 =8 
■ -<£ 2 >^ri=8 andr 2 =8 

• {i?>->n=8 and ^=8 

Action Density -> rj=0 and r 2 =8 
▲ Eneigy Density ri =0 and r 2 =8 
r\=Q and ^=8 

• (E 2 ) -> n=0 and r 2 =8 



Figure 13: Comparison between the results for (r\ = 0;r2 = 8) and (r\ = 8;r2 = 8) along segment 
gluon-antiquark in the L geometry. 
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Figure 14: String attraction and fusion, and string repulsion, respectively in type I and II supercon- 
ductors. 



4. A double vertex model for the scalar confining interaction 

Our lattice computation of flux tubes, string-like extended objects, suggest that the confining 
potential is scalar since the groundstate of a bosonic string is a scalar. This is also suggested by 
the spectroscopy of hadrons. In a pertubative QCD scenario, the hadron spectroscopy would be 
qualitatively similar to electronic spectra of the lighter atoms. It is remarkable that the Spin-Orbit 
potential (also called fine interaction in atomic physics) turns out to be suppressed in hadronic 
spectra because it is smaller than the Spin-Spin potential (also called hyperfme interaction). This 
constitutes an evidence of non-pertubative QCD. Another evidence of non-pertubative QCD is 
present in the angular and radial excitations of hadrons, which fit linear trajectories in Regge plots, 
and suggest a long range, probably linear, confining potential for the quarks. This led Henriques, 
Kellet and Moorhouse, Isgur and Karl, and others [7, 8] to develop a Quark Model where a short- 
range vector potential plus a long-range scalar potential partly cancel the Spin-Orbit interaction. 
The short range potential is inspired in the one gluon exchange, and the quark vertex is a Coulomb- 
like potential, with a vector coupling \j/y^ y/. The long range potential has scalar coupling xj/y/, and 
is a linear potential. 

On the other hand, the QCD Lagrangian is chiral invariant in the limit of vanishing quark 
masses, because the coupling of a gluon to a quark is vector like. This is crucial because spon- 
taneous chiral symmetry breaking is accepted to occur in low energy hadronic physics, for the 
light flavors u, d and s, where, m M ,mj«m 5 < Aqcd < M N /3. The techniques of current algebra 
lead to beautifully correct theorems, the PCAC (Partially Conserved Axial Current) theorems. The 
Quark Models are widely used as a simplification of QCD, convenient to study quark bound states 
and hadron scattering. Recently [6] we have shown that these beautiful PCAC theorems, like the 
Weinberg theorem for n — n scattering, are reproduced by quark models with spontaneous chiral 
symmetry breaking. 

We aim to couple a quark line in a Feynman diagram with a scalar string, using the vector 
gluon-quark coupling of QCD. We remark that a simple vertex does not solve this problem, there- 
fore we reformulate the standard coupling of the quark to the confining potential. The coupling 
needs at least a double vertex, similar to the vertices that couple a quark to a gluon ladder in mod- 
els of the pomeron. Our double vertex is depicted in Fig. 15. 

To get the Dirac coupling of each gluon to a light quark, we follow the coupling obtained in 
the heavy-light quark system, computed in the local coordinate gauge [9]. This results in a Dirac 
coupling a with pair of y° matrices, which is also compatible with the Coulomb gauge. 
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r r 




V(\p-q\) 

Figure 15: The coupling of a quark to a string with a double gluon vertex 

In the color sector, each sub-vertex couples with the Gell-Mann matrix X a /2. Moreover the 

string is also a colored object because it contains the flux of color electric field. In quark models 

the string usually couples with a X a /2 to the quark line, here it couples with the two X a /2 of the 

sub-vertices. For a scalar coupling, which is symmetric, we use the symmetric structure function 
d abc defined with? 

{X\X b } = d abc X c . (4.1) 




Figure 16: The self-energy term of the mass gap equation 

The dependence in the relative momentum must comply with the linear confinement which is 
derived from the string picture, 

-'"<»-«> ~'£ji^F (42) 

where <7 ~ 200MeV /Fm. is the string constant, and ^ is an algebraic color factor. 

The gluon propagators and the different sub- vertices result in a distribution in the loop mo- 
mentum k, see Fig. 15. Here different choices would be possible. For simplicity we assume that 
the same momentum flows in each of the two effective gluon lines. We also remark that the 
distribution in k is normalized to unity once the correct string tension is included in the relative 
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potential V(p — q). This amounts to consider that the momentum k distribution is a Dirac delta, 

p + q 



(2?r) 3 5 3 ( k 



(4.3) 



We finally compute the vertex, decomposing the Dirac propagator in the convenient particle 
and anti-particle propagators, computing the energy loop integral, and summing in color indices, 



% ff = x c (s k +c k k-y) 



(4.4) 



k= p+v 



where S k = m k j \/k 2 + m k 2 , C k = kj \/k 2 + m k 2 , and where m k is the constituent quark mass, to be 
determined in the next section. 



5. Results for the model of scalar confinement 
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Figure 17: Testing the convergence of the numerical method with the quark condensate (y/y/). 



We derive the mass gap equation projecting the Schwinger-Dyson equation with Dirac spinors, 

u s (p)y -\p)v s ,(p) -u s (p)Z(p)v s ip) = (5.1) 

where J?b is the free Dirac propagator, and where £(/?) is the self-energy depicted in Fig. 16. The 
self energy consits in a three loop Feynman diagram, including one loop in each double vertex and 
a third rainbow-like loop for the string exchange interaction. Nevertheless each double vertex is 
simple, see eq. (4.4), and we get for the self-energy term, 



u s (p)L(p)v s f(p) = J [(Cf -S 2 k )S q C p + 2S k C k S q S p k-p- 



-C q S p q-p-2S k C k C q C p k-q + +2ClC q S p k- qk- p]V (\p - q\) . (5.2) 
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The mass gap equation is a difficult non-linear integral equation, that does not converge with 
the usual methods [10, 11, 12]. Here we develop a method to solve it with a differential equation, 
using a convergence parameter X — ► 0. Our technique consists in starting with a large infrared 
cutoff X, where the integral term in eq. (5.2) is small. Then eq. (5.1) for the chiral angle (p p is 
essentially a differential equation which can be solved with the standard shooting method [12]. 
Next we decrease step by step the X parameter, using as an initial guess for the evaluation of the 
integral the (p p determined for the previous value of X. We test the convergence of the method 
computing the quark condensate (yn//), see Fig. 17. 




We solve the mass gap equation for the spontaneous breaking of chiral symmetry. In Fig. 
18 we compare the constituent quark mass m^, computed with our double vertex defined in eq. 
(4.4), with the mass computed with a simple Coulomb gauge [1 1] vertex X c y°. It turns out that the 
dynamical quark mass computed with the double vertex is closer to the one estimated is QCD 
sum rules, and this is a good point for the present work. 

In the chiral limit of a vanishing quark mass, the effective vertex (4.4) %ff—> X c k • y is pro- 
portional to y M and is therefore chiral invariant as it should be, whereas in the heavy quark limit 
%ff — > X c is simply a scalar vertex. The dynamical generation of a quark mass also generates 
a scalar coupling for light quarks. 

6. Conclusions 

We investigate the components of the chromoelectric and chromomagnetic fields around a 
static gluon-quark-antiquark [13]. We study the fields with two geometries, U and L shaped for 
the Wilson loops. This allows us to determine the shape of the flux-tubes with respect both to 
the energy and action densities. When the quark and the antiquark are superposed the results are 
consistent with the degenerate case of the two gluon glueball, and when the gluon and the antiquark 



14 



Schwinger-Dyson equations and the quark- antiquark static potential 



Pedro Bicudo 



are superposed the results are consistent with the static quark-antiquark case. The general shape of 
the flux-tubes observed here reinforce the conclusions of Bicudo et al. [2] and Cardoso et al. [3]. 
Another interesting result is that the absolute value of the chromoelectric field dominates over the 
absolute value of the chromomagnetic field. 

Moreover, in this talk we also present a model for the coupling of a quark to a scalar string [14]. 
Since the string is composed of bosonic gluons, it is expected that the groundstate of the string, 
dominating confinement, is a scalar. To model scalar coupling, double vector vertices are used, and 
the quark confining interaction has a single parameter, the string tension a c± 200MeV /Fm. The 
results are encouraging, our model not only enhances chiral symmetry breaking, but the matching 
from vector confinement at quark high momenta and scalar confinement at small quark momenta 
is beautifully achieved. A similar matching has been obtained by Alkofer, Fischer, Llanes-Estrada 
and Schwenzer [15]. We will try in the future to reproduce the whole hadron spectrum. 
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